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Developments and new applications of numerical stochastic perturbation 
theory 
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A review of new developments in numerical stochastic perturbation theory (NSPT) is presented. In particular, 
the status of the extension of the method to gauge fixed lattice QCD is reviewed and a first application to compact 
(scalar) QED is presented. Lacking still a general proof of the convergence of the underlying stochastic processes, 
a self-consistent method for testing the results is discussed. 
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1. NSPT for gauge fixed lattice QCD 

Numerical stochastic perturba- 

tion theory (NSPT) was successful in performing 
high loops gauge invariant computations in lat- 
tice QCD (for new, and potentially relevant, 
phenomenological implications see ||). An exten- 
sion to gauge fixed computations would of course 
Le highly welcome. First results in this direction 
were reported in St. Louis, but the matter turned 
out to be not settled down. In the following, as a 
first step, the basics from ||| are recalled. 

As it is well known, only gauge invariant quan- 
tities have a large time limit in standard stochas- 
tic quantization jj], while gauge-non-invariant 
quantities are affected by divergences. As usual, 
the situation is slightly different on the lattice due 
to the fact that the gauge degrees of freedom have 
been compactified, so that gauge-non-invariant 
quantities average to zero. Gauge fixed Langevin 
simulations rely on the old observation || that 
the standard evolution step 
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(77 is the gaussian noise) can be interlaced with a 
gauge transformation 
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the convergence of the process for gauge invariant 
quantities is unaffected, while the gauge trans- 
formation G can be chosen in order to enforce 
the gauge condition one is interested in: the 
stochastic process is then attracted towards the 
submanifold defined by the gauge condition and 
gauge-non-invariant quantities no longer average 
to zero. This recipe is in a sense the lattice imple- 
mentation of Zwanziger's stochastic gauge fixing 
scheme [||. 

NSPT rely on an expansion of the solution 
of Langevin equation in powers of the coupling, 
which actually decompactifies the fields (in a 
sense the fundamental fields of NSPT live in the 
Lie algebra), so that divergences in the gauge- 
non-invariant sector are back. Nevertheless one 
can adopt the stochastic gauge fixing scheme also 
in this context, with the caveat that also the 
gauge condition one wants to enforce has to be 
expanded as a series of conditions. To be defi- 
nite, as the field is expanded as = A^ + 
gA^jP +g 2 A^ + . . ., the form of Landau condition 
one needs to impose is 

d,A^ = (4) 
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for every order k (partial derivatives, as usual, are 
to be understood as finite difference operators). 
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The implementation presented in ||] goes as fol- 
lows. One considers the functional (norm) 

N[G] = X>(A°(n)Af (n)) = $> fc jV (fc) (5) 

where an expansion for N is induced by the ex- 
pansion of the field A^. If the gauge transforma- 
tion is chosen as 

G = e- aw w = Y J 9 k d^ (6) 

k 

extremum conditions for every are recov- 

ered that exactly enforce ([I]). First simulations 
with this recipe were presented in |jj and they 
seemed good up to order g 4 . As a matter of fact, 
if one lets the system evolve for a while, the signal 
for gauge-non-invariant quantities is not at all 
stable, but one recovers a collection of plateaux. 
While gauge fixing transformation is in charge to 
kill divergences due to "longitudinal" degrees of 
freedom, (||) lets some divergence still at work. 
This can be checked by inspecting the norms of 
the fields . While the norm of A u is mini- 
mized and stable through the process (as a mat- 
ter of fact, this norm is exactly N^ 2 \ which is 
the only ATw one is ensured to minimize) , higher 
orders suffer divergences in their norms. 

A first attempt has been done in order to 
enforce (|j) and keep every norm under control 
through the process. This can be obtained by 
adding new terms to «/ fe ' in (g) (the formulae are 
quite cumbersome and they will be written down 
elsewhere). The results are encouraging and can 
be summarized by inspecting fig. [j]. In the up- 
per half the norm of A^ is monitored through 
the evolution with and without these new terms 
((Q) is enforced in either case). The lower half 
shows the signal for order g 4 of the trace of the 
link, again with and without these new terms: 
the new gauge fixing procedure keeps the plateau 
stable (while the old procedure lets the system 
depart from it). 

2. How can you trust NSPT? 

Of course a stable signal does not mean at all 
that a result is correct. Actually there is an em- 
barrassing lack of a general proof of the conver- 
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Figure 1. Above: the norm of the field A$ with 
and without corrections to @. Below: the signal 
for order g 4 of the trace of the link in the same 
simulations. 



gence of the (quite cumbersome) stochastic pro- 
cesses underlying NSPT, so that in a sense one 
could worry about whether to trust a NSPT re- 
sult or not. A general strategy to check a result 
can be outlined as follows: check up by explicit 
computation the trivial order and then set up re- 
lations among different perturbative orders. To 
be definite, one can think of Schwinger-Dyson 
equations as obtained by 

- J JtyA(OMe-*M) = (7) 

By expanding this equation in the coupling one 
recovers relations among different orders. As it is 
well known, this is a standard approach to per- 
turbation theory. In this context one wants to 
look at these relations to test whether they holds 
once every expectation value involved has been 
obtained by means of the Langevin technique. 
Unfortunately for (lattice) QCD ([?]) reads 

- J DUV{0[U]e- Sw[u] e- SGF[u] k FP ) = (8) 
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where V is the group derivative and App the 
Fadeev-Popov determinant, which is tremendous 
to take into account without introducing ghosts. 

Nevertheless, from the standard formula for a 
QCD expectation value 

f DADcDc 0[A] e~ s ^ e~ s ^ A ^ (9) 

Z{g) J 

(where the action has been split in two terms S 
and S c , the latter being dependent on the ghosts), 
one can get a recursive relation for the coefficients 
of the expansion ~^2 k g k O^ . Namely, by deriving 
with respect to the coupling <?, one gets (O' means 
d/dg O) 

(k + l)0 {k+1 ^ = O'W - (OS") (fc) 



l-\-m—k 



(10) 

+ ]T ^' (0 o (m) - (os c 'Y k) 

l+tn—k 

One can obtain by the Langevin technique ev- 
erything except the last two contributions, which 
depend on the ghosts. For example, testing the 
order g 4 result for the trace of the link only re- 
quires to compute the ghost corrections to the 
gluon propagator. 

3. NSPT for compact QED 

A first implementation of the Schwinger-Dyson 
test program has been done within the context of 
a new application of NSPT, namely to compact 
QED. The final goal of this project is a collab- 
oration with the Tor Vergata group, aiming at 
merging NSPT with the bermions approach 
to push one loop further the computation of elec- 
tron's g — 2. An important point to notice is that 
the extrapolation in the number of flavors is per- 
turbatively well defined. Being QED not asymp- 
totically free, a major problem to deal with is an 
explicit inclusion of counterterms. 

In order to gain experience, we started simu- 
lating scalar QED defined by the action 

S= ^l](l-cos^)-^0*(n)(C/ M (ri)^(n+A) 

P n)i 

+Ul(n - AM" - /})) + Y( M2 + 8)f(«)^) 



<^ t „ being defined by the plaquette U uv = 

We have till now performed computations up 
to order a 3 in the gauge invariant sector, every 
result having been tested by Schwinger-Dyson 
equations. In many instances the overall picture 
is the same as for QCD; in particular, some form 
of stochastic gauge fixing is needed also in the 
gauge invariant sector in order to keep fluctua- 
tions under control. 

4. More to come . . . 

Another application which is expected to come 
soon (some preliminary work has just started af- 
ter the conference) is a high order perturbative 
expansion of the QCD running coupling asF as 
defined in the SF scheme ^ . One is interested in 
pushing further the connection between different 
schemes and in evaluating the high order behavior 
itself of this definition of the coupling (in Q the 
coupling a p i defined via the plaquette has been 
found to be affected by a renormalon). 
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